Long quantum transitions due to unstable semiclassical 

dynamics 

O ■ D.G. Levkov: A.G. Panint 

o ■ 

^ , Institute for Nuclear Research of the Russian Academy of Sciences, 

, 60th October Anniversary Prospect 7a, Moscow 117312, Russia 

o : 
^ . 

Abstract 

\ Quantum transitions are described semiclassically as motions of systems along 

I ■ (complex) trajectories. We consider the cases when the semiclassical trajectories are 

^ , unstable and find that durations of the corresponding transitions are large. In ad- 

\ dition, we show that the probability distributions over transition times have unusual 

qh' asymmetric form in cases of unstable trajectories. We investigate in detail three types 
'— of processes related to unstable semiclassical dynamics. First, we analyze recently 

\ proposed mechanism of multidimensional tunneling where transitions proceed by for- 

J> ' mation and subsequent decay of classically unstable "states." The second class of 

■ processes includes one-dimensional activation transitions due to energy dispersion. In 
Q>^ \ this case the semiclassical transition-time distributions have universal form. Third, 

' we investigate long-time asymptotics of transition-time distributions in the case of 

ff^ ' over-barrier wave packet transmissions. We show that behavior of these asymptotics 

, is controlled by unstable semiclassical trajectories which linger near the barrier top. 

o : 

: 1 Introduction 

X: 

■ Semiclassical method unveils fascinating relation between the properties of quantum transi- 
tions and the character of classical dynamics of the system. This relation is valid, in partic- 
ular, for tunneling processes [l] which do not exist at the classical level. Indeed, tunnelings 
are described semiclassically by complex trajectories [2] evolving classically in complex phase 
space. Observables characterizing tunneling processes, such as probabiUties [2] or splittings 
of energy levels [3], are functionals of complex trajectories. Clearly, features of (complex) 
classical dynamics are imprinted in these observables. 

The link between tunneling processes and classical dynamics is supported by theoretical 
studies of multidimensional tunneling. It was found that behavior of tunnel splittings of 
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Figure 1: (Color online) Experimental setup for measuring the traversal-time distribution 
p(r). Only the "transition" coordinate x is shown. Macroscopic detector is represented by 
the hatched region x > X2- 

energy levels is qualitatively different for systems with integrable [3], near-integrable [5], 
mixed [6] and chaotic [7] dynamics. Some of the related phenomena have already been 
observed experimentally [8j. 

In this paper we study impact of non-trivial semiclassical dynamics on temporal char- 
acteristics of classically forbidden transitions. We are particularly interested in cases when 
complex trajectories describing transitions are unstable. Such trajectories do not interpolate 
directly between in- and out- regions of the process, but rather get attracted to a certain 
unstable orbit lying on the boundary between the regions. It takes large or even infinite 
time for a trajectory of this kind to arrive to the out-region. In general one does not expect 
the time scale of the semiclassical motion to coincide with the (properly defined) duration 
of quantum transition. Nevertheless, we advocate qualitative relation [9] between these 
quantities: durations of quantum processes are large when the corresponding semiclassical 
trajectories are unstable. 

Presently there is no generally accepted definition of quantum traversal time pO]. During 
last decades a number of definitions have been suggested in literature [HI [12], [131 [El [15] . 
Below we follow the approach of Refs. IT7[ [18] where traversal time r is treated as 
a random variable with probability distribution p(t). Experiment for measuring p(r) is 
depicted schematically in Fig. [TJ The horizontal axis in this figure represents "transition" 
coordinate x which runs between in- and out- regions of the process and meets potential 
barrier f/(x) on the way. Quantum particle is emitted at t = 0, x ~ Xi in a definite initial 
state ^I/i(x) = \I'(0,x). As time goes on, the particle moves towards the potential barrier 
and eventually traverses it. At t = r the particle is caught by the asymptotically distant 
macroscopic detector which occupies the region x > X2 and registers the time of arrival r. 
The probability of registration is 

Vr= j d^j\m{T,^j)\\ (1) 

X>X2 

where the integration variable is denoted by xj for later convenience. Note that the total 
transmission probability Voo is obtained from Vt in the limit r +oo. We also remark that 
the points xi and X2 are assumed to be far enough from the potential barrier, so that the 
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particle moves freely at x < xi, x > X2. One introduces the normalized probability density 



By construction, p(r) ■ dr is the increment of the registration probability, which is conditional 
with respect to the transition event. Following Refs. ^T6\ \T7\ ITS] , we associate p(r) with 
probability distribution over traversal times r. The average time of passing is then computed 

as 



and time dispersion is defined accordingly. 

It is worth noting that several assumptions about detector registering the particles have 
been made in Eq. ([T]). First, we adopted approximation of instantaneous measurements 
disregarding time resolution of the detector. Second, we assumed that the states of particles 
are not distorted prior to the measurement event. Third, our probability formula ([1]) is 
applicable only in the case of small registration probabilities: true probability of registering 
the particle at t = r is conditional with respect to previous non-detection. For classically 
forbidden transitions considered below the probability of non-detection is close to 1 and 
Eq. (pP) is applicable. In other cases expression for Vr should be correctecH [TP] . 

Equipped with the definition ([3]), we study three examples of quantum processes related 
to unstable semiclassical trajectories. In all cases our results confirm the qualitative conclu- 
sion [9] about large transition times. Namely, we show that (r) ~ | logg'^\, o"^ > 0{g^), where 
g'^ is the semiclassical parameteiE|. In the limit of vanishingly small g^ when the system be- 
comes increasingly "more classical" the average duration of transition grows to infinity, while 
0"^ stays finite or grows. This behavior is drastically different from the respective scalings in 
the case of stable trajectories where (r) reaches definite limit and vanishes as g^ 0. 

We consider the following classically forbidden processes. First and most importantly, we 
study tunneling transitions in non-linear multidimensional systems at high energies. Insta- 
bility of semiclassical dynamics in this case was demonstrated in Refs. [201 EI], see Refs. [22j 
for earlier works. We have already sketched the behavior of the corresponding complex tra- 
jectories: they approach certain unstable orbit, remain close to it for an infinite time interval 
and then slide away into the out-region. In systems with two degrees of freedom the media- 
tor unstable orbit describes periodic oscillations near the saddle point of the potential. We 
call it sphaleroT^ or simply unstable periodic orbit. We also attribute instability of respective 

^There are two important simplifications in our setup. First, we prepare definite in-states of the process. 
Second, we consider the case of asymptotically distant detector, which saves us from the controversial issue 
of separating transmitted and reflected parts of wave function. 

^One can simply consider detectors registering particles with probabilities i^Vr, where v is a. small number. 
In this case Eq. ^ is valid, since the probabilities of previous detections are small. 

•^This is a dimensionless quantity proportional to h which measures the relative size of quantum fluctua- 
tions. Semiclassical expressions become accurate as — s- 0. 

''This term is taken from field theory [23]. We find it more convenient than exact but lengthy expression 
"normally hyperbolic invariant manifold" (NHIM) |24) . 
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trajectories to appearance of new mechanism of sphaleron- driven tunneling (or tunneling 
driven by unstable periodic orbit). Recent studies [25l[26l[27j revealed two experimental sig- 
natures of this mechanism: the probabilities of related processes are suppressed by additional 
power-law factor and distributions over out-state quantum numbers are anomalously 
wide pn 1^ Below we show that the distribution p{t) and corresponding values of (r) 
and 0"^ can be used for experimental identification of the sphaleron-driven mechanism of 
tunneling. 

Second, we investigate wave packet transmissions through one-dimensional potential bar- 
rier. We consider the case when the average energy of wave packet is lower than the height 
of the barrier but sufficiently close to it. Then transmission proceeds by activation caused 
by energy uncertainty in the initial state. The probability of such transmission is expo- 
nentially suppressed, while the corresponding semiclassical trajectory approaches the top of 
the barrier and stays there for an infinite time. This trajectory is unstable; its behavior is, 
again, of the type described abov^. Our calculations show that the forms of the semiclas- 
sical traversal-time distributions are universal in one dimension: all information about the 
potential barrier and initial wave packet is encoded in two parameters of the distributions, 
which can be absorbed by shifting and rescaling r. Universal forms of p(r) lead to universal 
qualitative dependencies (r) ~ | logg^\, cr^ ~ 0{g'^) in the case of one-dimensional activation 
transitions. 

Third and finally, we study the probabilities of large time delays r in one-dimensional 
over-barrier transmissions. In this case the average energy of particle exceeds the height 
of the potential barrier. The particle is registered at time r after transmission at a given 
distance behind the barrier; we calculate semiclassically the asymptotics r +00 of the 
registration probability. Semiclassical trajectories entering this calculation linger near the 
barrier top. The probability at large r has universal form, again. Due to delays of the 
respective trajectories this probability is larger than one would expect naively. 

We illustrate our findings by performing explicit calculations in systems with one and 
two degrees of freedom. In these systems we compare the semiclassical results for p(r) with 
results of explicit quantum-mechanical calculations. We find agreement at small g^. 

The paper is organized as follows. We start in Sec. [2] with properties of traversal-time 
definition Q. Then, in Sec. Owe discuss the mechanism of sphaleron-driven tunneling and 
one-dimensional activation processes. General semiclassical expression for p(r) is introduced 
in Sec. HI Results for the durations of quantum processes are presented in Sec. |5l we consider 
one-dimensional activation transitions in Sec. 15. ![ investigate sphaleron-driven tunneling in 
Sec. 15.21 and study long-time asymptotics of "over-barrier" traversal-time distributions in 
Sec. 15.31 We summarize in Sec. O Details of semiclassical and exact quantum-mechanical 
calculations are presented in Appendices. 

^The constituent unstable orbit in this case is static classical solution "sitting" on top of the barrier. 
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2 Properties of traversal— time definition 



The controversial issue of quantum traversal time has attracted considerable attention dur- 
ing last decades (TU]. The reason for the controversy seems to be hidden within the quantum 
theory itself which does not offer clear definition or unique method of computing the dura- 
tions of quantum transitions. At the moment there exists a number of physically different 
definitions of quantum traversal time, see e.g. Refs. [HI [121 [13 fl^ [15] . 

Expression ([3]) belongs to a wide class of definitions involving delay of the transmitted 
wave packet with respect to the incident one. In this approach the delay is computed using 
some feature of wave packet, maximum [TTj or centroid [2S]. In our case the probability 
conservation law implies that the distribution ([2]) is proportional to the total current through 
the surface x = X2 at t = r. Thus, p(r) reflects the form of the transmitted wave packet, 
and mean time of passing Q fits into the above class of wave packet -related definitions. We 
remark that p(r) ■ dr can be regarded as a probability of registering the particle at time r 
by the point-like detector situated at x = X2. 

The properties of the definition ([3]) are similar to those of other definitions from the 
same class. First, the setup used in the Introduction resembles recent experiments with 
photons [291 [30], cind one can hope to measure the distribution p(r) experimentally. Note, 
however, that experimental verification of our results, if possible, requires substantial mod- 
ification of existing experiments. Indeed, recent measurements of traversal times have been 
performed in quasi-stationary regime where wave functions of particles are substantially 
wider than potential barriers. In this case the interval Q coincides [281 [H] with the phase 
time of Bohm and Wigner [11] which is reproduced in actual measurements [29]. On the other 
hand, our results are obtained in the "semiclassical" case where the coordinate uncertainty 
is parametrically smaller than the barrier width. 

Second, the interval passes [IB] self-consistency check of Refs. [lOj: times of trans- 
mission through the barrier and reflection from it have sense of conditional averages over 
two mutually exclusive events. 

Third, our traversal time ([3|) suffers from Hartman effect [31] which leads to superluminal 
velocities of under-barrier motions. This effect was observed experimentally [SO]- It is 
explained as follows. One notices [IHl [22] that the out-state of the tunneling process is 
formed by the waves constituting forward tail of the initial wave packet. Thus, the out-state 
is not related casually to the central part of the in-state wave function. On the other hand, 
only the central parts of in- and out- states contribute substantially into (r). Since casual 
connection between these parts is absent, apparent velocity {x2—Xi)/ (r) can be superluminal. 
We stress that the definition is sufficient for the purposes of this paper: we regard (r) as 
an experimentally measurable [30j quantity characterizing quantum transitions. 

Finally, let us remark that quantum traversal time can be naturally defined without 
any reference to the delay between transmitted and incident wave packets: one can exploit 
stopwatches attached to the particle [T3] or observe response of total transmission probability 
to the periodic modulation of the potential [14j. In these approaches the properties of 
traversal times are essentially different from those listed above. Although we believe that the 
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qualitative conclusions of this paper should be valid for any reasonable setup, our quantitative 
results are not applicable for other traversal-time definitions. 

3 Three processes 

Let us describe three types of processes related to unstable semiclassical dynamics. 

We start with one-dimensional activation transitions of quantum particles through po- 
tential barriers U{x). In what follows we consider the setup in Fig. [H The initial state 
is chosen to be a Gaussian wave packet with central position xi, average momentum 
Po and momentum dispersion a^. Transmissions proceed between the in- and out- regions, 
X —oo and x —>■ +00 respectively. Points xi and X2 in Fig. [1] belong to these regions. 

We exploit the semiclassical approximation which is justified if 

^2 ^ ^ 

Here Vq and /q are height and width of the potential barrier and m is the particle mass. 
Dimensionless parameter g"^ introduced in Eq. (jlj) measures accuracy of the semiclassical 
approximation; one can adjust the value of this parameter by changing U{x). Below we 
adopt dimensionless units m = h = Vq/II = 1. In these units the height Vq = l/g"^ and 
width Iq = 1/g of the "semiclassical" barrier are parametrically large. 

We use the following scalings of the in-state parameters with g"^. First, we suppose that 
the typical energy of the process is large, Po ~ ^0 ~ ^/d^- Second, momentum dispersion 
is chosen to be of order one, 0"^ ~ 0{g^). Finally, the points Xi and X2 are taken to be 
sufficiently far from the potential barrier, so that motion is free at x < a;i, x > X2. This 
implies xi, X2 ~ l/g. 

We have already mentioned in the Introduction that one-dimensional transmissions pro- 
ceed by activation at some values of in-state parameters. Let us explain the origin of the 
effect by decomposing the initial state \Efj in the basis of plane waves with fixed momenta 
p. We consider the case po < Pc,2 = -\/2Vo when the initial wave packet is mostly formed 
by "under-barrier" modes with p ~ po < V2Vo. Clearly, these modes are exponentially 
suppressed after transmission. On the other hand, "over-barrier" modes with p > v^2T^ 
constitute exponentially small fraction of the initial state, but overcome the barrier with 
probability of order one. One sees that, depending on the suppressions, the probability 
Voo can be dominated by "under-barrier" or "over-barrier" modes. We call the respective 
transitions by tunneling and activation. Note that the value of Voo is exponentially sup- 
pressed in both cases. Note also that the activation processes are related to the possibility 
of overcoming the barrier due to finite momentum dispersion in the in-state. 

In Sec. 15.11 we will argue that there always exists some critical value of momentum pc,i 
such that transitions proceed by activation in the region pc,i < Po < Pc,2- We will also 
demonstrate explicitly that activation transitions are described by unstable semiclassical 
trajectories which linger near the barrier top. 
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Figure 2: (Color online) Potential (not to scale). On top of the potential surface: real part 
of complex trajectory (thin red line) describing sphaleron-driven tunneling at gpo = 1.55, 
o"p = 0, g'^Ey = 0.05 and the sphaleron orbit (thick black line). The saddle point of the 
potential is shown with thick black dot. 

The second type of processes is related to the multidimensional mechanism of sphaleron- 
driven tunneling. For definiteness we explain this mechanism in the model of Refs. [331 [201 [2S] 
where two-dimensional quantum particle of mass m moves in harmonic potential valley of 
frequency ujq. The valley is intersected at an angle by the potential barrier. The potential 
of the model is 

U{x,y) = ^y^ + Voe-^^^y^'/^^o, (5) 

where Vq and /q are height and width of the barrier, respectively. Function U (x, y) is plotted 
in Fig. [1 

Semiclassical approximation in the model (jSj) is justified under conditions (jlj) and 

u = uJo lo^/m/Vo ~ 0{g°) . (6) 

In dimensionless units introduced above the potential (jSD contains only two parameters, g^ 
and u, 

U{x,y)=uV/2 + \e-^'^^+y^'/' . (7) 
9^ 

We use uj = 1/2 for the waveguide frequency. 

In what follows we study tunneling transitions between the asymptotic regions x ^ —oo 
and X — * +00 (in- and out- regions respectively). In-states of transitions are fixed 
completely by taking \I'j(x,?/) = ipx{x) ■ ipy{y), where ipxi^) is a Gaussian wave packet with 
central position xi, average momentum pq and momentum dispersion cXp, and ipy{y) is an 
oscillator eigenfunction with energy Ey. We propagate semi classically the in-state into the 
out-region and compute the probability Voo and traversal-time distribution p(r). 
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In the semiclassical regime the typical energies of motion in y direction are parametrically 
large, Ey ~ 1/ g"^. The scalings of other parameters are chosen in the same way as in the above 
case of one-dimensional transitions. We characterize transitions with the combinations gpQ, 
g'^Ey, (Tp, gxi, gx2 which stay finite as (7^ — 0. 

Note that transmissions between the asymptotic regions of the potential ([7]) remain clas- 
sically forbidden even when the total energy pQ/2-|-ii^y exceeds the minimum height Vq = 1/g'^ 
of the potential barrier. Indeed, classical trajectories with energies slightly higher than Vq 
can overcome the barrier only by passing in the immediate vicinity of the saddle point of the 
potential. General trajectory misses this vicinity, and the corresponding transition remains 
classically forbidden. The probability of such transition is exponentially suppressed on gen- 
eral grounds. Tunneling transitions at energies exceeding the height of the potential barrier 
are known as dynamical tunneling [21 [M] . 

The mechanism of sphaleron-driven tunneling is based [201 EI] on behavior of complex 
trajectories in the regime of dynamical tunneling. One of such trajectories is depicted in 
Fig. [2] (thin red line). It overcomes the barrier in two stages, by approaching the unstable 
periodic orbit (sphaleron) and departing from it. Note that the sphaleron is real; it is plot- 
ted in Fig. [2] by thick black line. This orbit is unstable by construction, since it describes 
oscillations around the saddle point x = y = of the potential. The semiclassical tra- 
jectories in the new tunneling mechanism remain close to the sphaleron for an infinite time 
interval; imaginary parts of these trajectories become negligibly small after several sphaleron 
oscillations. 

In general one finds that the sphaleron-driven mechanism is relevant in some region of 
in-state parameters, pc,i < Po < Pc,2, where pc,i, Pc,2 are critical momenta depending on Ey 
and (Tp. Regions of low and high momenta, po < pc,i and po > Vc,2 respectively, correspond 
to direct tunneling through the barrier and classical over-barrier transitions. 

We remark that multidimensional processes of sphaleron-driven tunneling are essentially 
different from the one-dimensional activation transitions considered above: instability of 
trajectories in the sphaleron-driven case is not caused by finite momentum dispersion. 

The third effect considered in this paper is related to the possibility of large time delays in 
over-barrier wave packet transmissions. For simplicity we return to one-dimensional setup 
where Gaussian wave packets traverse the potential barrier U{x). Now, we are interested in 
the case of high initial momenta, po > -\/2Vo, when the central part of wave packet overcomes 
the barrier classically. We study the tail of p(r) corresponding to large traversal times r. In 
Sec. 15. 31 we will argue that this tail can be readily found by considering unstable semiclassical 
trajectories which linger near the barrier top. 

From the physical viewpoint the origin of the effect is explained, again, by decomposing 
the initial state into a collection of plane waves with fixed asymptotic momenta p. In general, 
waves with p ^ a/2Vo are populated, although the respective probabihty is exponentially 
small. These waves get delayed near the top of the barrier and create long tail of p(t). 
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4 Semiclassical method 



Here we describe the standard method of complex trajectories [2] and introduce convenient 
semiclassical expression for p(t). Although the semiclassical techniques of this Section are 
pretty general, we concentrate on the one- and two-dimensional models of Sec. [31 

For a start, let us briefly explain the origin of complex trajectories [2] leaving details 
of semiclassical computations to Appendix |Al Consider evaluation of total transmission 
probability Poo- One starts with the path integral for the final state, 

e^^M ^ (8) 

where the paths x(t) are defined in the time interval t G [0; r]: they start at x(0) = Xj and 
arrive to x(r) = x/. S[x\ is the action in this time interval. Since we are interested in Voo, 
the value of r will be sent to infinity in the end of calculation. 

At small the integral ([H]) can be evaluated by the saddle-point method. Complex 
trajectory x(t) describing tunneling process is the saddle-point path of integration. In 
Appendix |X] we obtain saddle-point equations for this trajectory. Namely, x(t) should 
satisfy the classical equations of motion. Conditions at t = come from the saddle-point 
integration over Xj, 

Xi - 2ial{xi - xi) = po , Ey = {{/i +uj'^yi)/2 , (9) 

where the subscript i denotes quantities computed at t = 0. [In the one-dimensional model of 
Sec. [3] only the first condition should be imposed.] One sees that the second of Eqs. fixes 
initial energy of y-oscillator. Interpretation of the first equation is clear if the trajectory 
x(i:) is real. In this case the real and imaginary parts of the first of Eqs. (Q constitute 
two independent conditions specifying the initial particle position Xi = xi and momentum 
Xi = Pq. In the complex case the two parts are combined with the coefficient proportional to 
the momentum dispersion a^. Final boundary conditions for x(t) arise when one substitutes 
\E'(r, xj) into the probability expression ([T]) and evaluates the saddle-point integral over xj. 
One obtains conditions 



^1 



r, X; 



rfXi \E'i,('x,; 



[dx] 



Imx(t), Imx(t) ^ as t +oo (10) 

which mean that evolution in the out-region is real. Note that the limit r —>■ +oo is already 
taken in Eq. ffTOj) : at finite r the integral over x/ cannot be evaluated by the saddle-point 
method. 

After finding the trajectory x(t) from the classical equations of motion and boundary 
conditions flTUl) . one computes the transmission probability by the semiclassical formulau 

Poo = A-e-^/^% (11) 
^We assume that only one trajectory gives substantial contribution into the probability. 
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where the suppression exponent F and prefactor A are functionals of the saddle-point tra- 
jectory x(t). Exphcit expressions for these functionals are derived in Appendices lAl and [C| 
see Eqs. (IST]), and (|63|). 

Now, let us evaluate semiclassically the distribution p(t). Following straightforward 
method [S], one should repeat the above considerations in the case of finite durations of 
transitions r and apply Eqs. ([I]), ([2]). It is simpler, however, to use indirect technique 
introduced in Refs. [26] . 

With the purpose of explaining the logic of Refs. |26], we start from purely classical case 
when the traversal time r is defined as the time of motion in the region x < X2- [Recall that 
real classical trajectories start at a; = Xi, t = 0.] One introduces the functional 



Tint [x] 



dt 6{x2 — x{t)) , 



(12) 



which returns traversal time for any classical path x(t) G M. We call Tint interaction time. 
To fix the value of Tint, one introduces Lagrange multiplier e and changes the classical action, 
^[x] S[yi] + e(Tj„([x] — r). Trajectories Xf(t) extremizing the new action correspond to 
fixed durations r = Tj„t[xj] of transitions. Note that the Lagrange method enables one to 
replace explicit fixation of traversal time r with modification of classical equations of motion. 

It is not trivial to introduce Lagrange method in quantum case. Indeed, there is no 
a priori reason to associate stochastic variable r with values taken by the functional Tint- 
Still, we find that there exists close connection between these quantities. The technical trick 
revealing this connection is explained as follows. One inserts into the integrand of Eq. ([H]) 
the unity factor 



cir'5(Ti„t[x] - t' 



dr' 

J too 



ide 
2V 



,E(r'-Ti„4x])/s2 



(13) 



Since the paths x(t) in Eq. (jH]) are defined in the interval t G [0; r], we temporarily restrict in- 
tegration in Eq. ( |T2l) to this interval. Now, Tint [x] takes values between and r; we explicitly 
used this fact in the first equality of Eq. (fT3l) . In the second equality we performed Fourier 
transformation and introduced "Lagrange multiplier" e. The path integral for \E'(r, xj) takes 
the form, 



v&Cr, X 



f) 



dr' 



ide 
2V 



dyii \E',;(xi) / [dx\ 



,iS[x]-£(T,„t[x]-T')/32 



(14) 



Recall that the paths x(t) in this integral are defined in the interval t G [0; r]. Expression 
is remarkable in three respects. First, its integrand depends on x(t) via the new action 
5'[x] + ze(Tj„t[x] — t')/(7^ which is similar to the action arising in the classical Lagrange 
method. This is natural, since by inserting the (5-function, Eq. f[T^ . we introduced explicit 
constraint Tj„f[x] = t' in the path integral. Second, apart from the modification of the 
action the integral in brackets in Eq. (fT^ has absolutely the same form as Eq. ([8]). Thus, it 
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is evaluated in the same way as in Eq. ([8]), by finding the modified complex trajectory Xe(t)El 
By construction, Xe(t) extremizes the modified action and satisfies the initial conditions IQ. 

The third favorable property of Eq. (fT^ is an explicit dependence of the integration limit 
on T. Let us schematically explain how this feature simplifies evaluation of the distribution 
p(r). Recall that p(r) is proportional to the derivative of Vr with respect to r, see Eq. 
This derivative is evaluated easily if we represent Vr in the form Vr = dr^ . . . , where r+ 
is an integration variable, and the integrand does not depend on r. Then, p(r) is simply 
proportional to the integrand of Vr taken at r+ = r. 

We derive the desired integral for Vr in Appendix [Bl We start from Eq. (HM and perform 
all saddle-point integrations except for the integral over r' which is kept intact until the last 
moment. Then, we substitute \l/(r, xj-) into the probability expression ([T]) and integrate over 
xj. In this way we obtain Vr in the integral form, where r explicitly enters the integration 
limits. At this stage the integrand of Vr is a functional of the modified trajectory X£(t) which 
is defined in the interval t G [0; r]. Using the modified equations of motion, we extend Xe(t) 
to t G [0; +oo] and show that motion at t > r does not contribute into the integrand of 
Vr- In this way we show that the integrand does not depend on r. Thus, p(r) is precisely 
the integrand of Vr at the end-point of integration; it is conveniently computed using the 
extended trajectory Xe(t), t G [0; +oo]. 

The result of the derivation sketched above (see Appendix [B] for details) is summarized 
as follows. One starts by finding the modified trajectory Xe(t) which satisfies equations of 
motion following from the modified action 

5,[x] = 5[x] + ^(T,„ax]-r) (15) 
9 

and old boundary conditions, Eqs. ([9]), ffTOj) . Importantly, Xe(t) is extended beyond the 
interval t G [0; r], and the final boundary conditions ffTOj) are imposed aX t —>■ +oo. Then, 
we determine the saddle-point value of e. One can show that e is real, so that the modification 
term in Eq. f[T^ is imaginary. This feature is specific to the case of classically forbidden 
transitions. Besides, e is related to the duration r by non-linear equation 

ReTint[yie] = r , (16) 

which is similar to the respective constraint in the classically allowed case. We remark in 
passing that e = corresponds to the original complex trajectory x(t) which saturates the 
probability Voo, see Eq. (fT5|) . 

After finding the modified trajectory Xe(t) and solving Eq. (fT6|) one computes the traversal- 
time distribution, 

p{r)=^^^-Ae-^^/^\ (17) 

^Since Xe(i) G C, the functional Tint should be continued into the complex domain. This is achieved by 
approximating 9(x) in Eq. (I12|) by a smooth function 6a{x); the width a of 6a{x) is sent to zero in the final 
semiclassical expressions. Note that the value of Tint is complex on complex paths. 
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where M is the normahzation factor. Suppression exponent and prefactor in Eq. ([1] 
are given by the same expressions as in Eq. (fTTj) . but with the modified action S'^ and modified 
trajectory Xe(t). 

Note that the factor A^ - ex.Y>{~Ft/ q'^) entering Eq. fll7p has meaning of "naive" transition 
probabihty computed for the trajectory with action (IT^ . The factor ^J—de/dr is not so 
trivial, however. From the technical viewpoint it arises due to fixation of Tint in the path 
integral ([8]). 

In what follows we exploit the relation which is derived in Appendix [Bl 

2e = -dFjdr . (18) 

Equation f|T8l) shows that the values of e and r are related by Legendre transformation. This 
is what one expects, since e is interpreted as a Lagrange multiplier conjugate to r. 

It is worth noting that the method of "Lagrange multipliers" introduced above resembles 
the approach of Refs. Note, however, one important difference: we define traversal- 

time distribution p(r) by Eq. ([2]) and then introduce the functional Tj„t[x] as a convenient 
technical tool for the semiclassical calculation of this distribution. The authors of Refs. [15] 
go the other way around. They start with the functional similar to Eq. (fT2l) and then define 
traversal time using path integral technique. 

Let us mention one important restriction of the semiclassical method of this Section. 
While deriving Eq. f|T71) we assumed that only one complex trajectory gives substantial 
contribution into p(t). This assumption is certainly correct for one-dimensional conservative 
systems; explicit calculations show that it also hold^S in the two-dimensional model ([7]). In 
other multidimensional systems [9j there might exist several (or even an infinite number of) 
trajectories giving comparable contributions into the traversal-time distribution. In this case 
Eq. ( IT7|) should be modified: contributions of different trajectories and respective interference 
terms should be accounted for. 

We finish this Section by considering the standard tunneling mechanism where the origi- 
nal complex trajectory x(t) is stable. We are going to show that in this case p(r) is symmetric 
and sharply peaked around the central value r = (r). From the Legendre transformation 0181) 
one sees that the exponent in Eq. f|T71) has minimum at e(r) = 0. This point corresponds 
to the original trajectory x(t). Since x(t) is stable, the respective value of r = ReTj^^x] 
is finite. One expands the exponent around the point e(r) = and obtains the Gaussian 
distribution, 

p.(r) = ^=-e-(-<-»V2.?^ (19) 



(Tt-v 27r 

where we restored the normalization factor and added subscript s (for "stable"). The central 



^In a wide interval r e At centered around (t) . Typically, the width of this interval is of order 40 and 
becomes somewhat smaller in the vicinity of the critical point pa ~ Pc,i- In what follows we work at r € At 
and use Eq. (fTT]) . 
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value (r) and dispersion of this distribution are given by the expressions 



(r) = ReT,„i[x]~0((?°) 



g^dr 
2 de 



(20) 



which come from the Taylor series expansion of the semiclassical exponent and Eqs. f[T6|) . 
f|T8|) . We remark that Eq. f|T9|) is applicable at |t — (r)| ~ where p(r) ~ 1. Outside of 
this region the value of p(r) is negligibly small. Note also that the average duration (r) in 
Eqs. f l2U]) is equal to the time spent by the original complex trajectory at x < X2. 

The above considerations are not valid if the original trajectory x(t) is unstable. In that 
case the point e = corresponds to r = +oo and the semiclassical exponent in Eq. (ITTl) 
remains almost constant as r — +oo. Thus, p(r) is essentially asymmetric: at large r it 
falls off slowly, only due to prefactor. In the next Section we will show that the corresponding 
values of (r) and cr^ are large. 



5 Results 

In this Section we compute distributions p(r) for three quantum processes related to unstable 
semiclassical dynamics. 

5.1 One dimension 

We begin by considering one-dimensional activation transitions introduced in Sec. [31 For 
definiteness we suppose that the potential U{x) reaches its maximum at a; = 0. We denote 
by a;^ = — f/"(0) the curvature of the potential at the maximum. Note that a;_ stays finite 
as (yf^ — > 0. We will also use the exemplary potential 

^^^^ g'^di^igx) ^ ^ 

for checking and illustrating the semiclassical technique. 

Let us show that the semiclassical dynamics is unstable in the case of activation tran- 
sitions. We use the logic of Sec. [3] and decompose the initial state in the basis of plane 
waves with fixed momenta p. Denoting by Vp the transmission coefficients of individual 
waves, we write, 

Poo = / dpVp- \^iip)\' = dpVp- j= e-(^'-^'«)'/2<x^ , (22) 

Jo Jo crW 27r 



where interference terms are absent due to conservation of energy Ep = p'^/2. Recall that the 
value of po is smaller than pc,2 = V2Vo. On the other hand, the integration momentum p in 
Eq. (!22l) can be higher or lower than y/2Vo representing the "over-barrier" or "under-barrier" 
modes, respectively. We consider these cases separately. Transmissions at p > \J2Vq are 
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described by classical over-barrier trajectories. Thus, Vp = 1 and the integrand in Eq. (l22l) 
decreases with p at p > y/2Vo due to the in-state contribution. 

At p < V2Vo the situation is different, since the coefficients Vp are exponentially sup- 
pressed. The semiclassical expression for these coefficients is well-known, Vp = exp[—S{Ep)] 
where S{Ep) is Euclidean action computed for one period of Euclidean oscillations in the 
upside-down potential. The semiclassical trajectory x{t) performing oscillations has energy 
Ep = Let us visualize how x{t) moves between the in- and out- asymptotic regions. 

In the in-region the trajectory evolves in real time until it stops at the left turning point 
of the potential. The subsequent semiclassical evolution proceeds in Euclidean time where 
x{t) performs precisely one half-period oscillation in the upside-down potential. At the 
end of oscillation the trajectory arrives at the right turning point. After that x{t) moves 
in real time, again, and finally ends up in the out-region. Only the Euclidean part of x{t) 
contributes into Vp. 

One sees that for under-barrier modes both factors in the integrand of Eq. ( !22l) are 
exponentially sensitive to the momentum p. To understand the behavior of the integrand, 
we consider its logarithmic derivative, 

^In [Vp\nP)\'] = pTiEp) - ^ , (23) 

where T{E) = —dS/dE is the period of Euclidean oscillations. Now, we take the limit 
p V^Vq in Eq. fl23|) . In this limit Ep approaches the height of the potential barrier, 
and Euclidean oscillations get restrained to a small region near the barrier top. Thus, 
T 2tt/uj_. Using this property, one finds that the derivative (1251) is positive at p = \/2Vq 
if Pq > Pel, where 

Pel = V^(l-27ra2/cu„) (24) 

is the critical momentum. Note that in the limit Ep — > Vq the semiclassical trajectory x{t) 
becomes unstable: it starts in the in-region and ends up on top of the potential barrier as 

t +00. 

Now, one recalls that the integrand in Eq. fl22|) decreases with p at p > ^/2Vq. Thus, it 
has sharp maximum atjo = ^/2Vq whenever po lies between pc,i and pc,2 = \/2Vo. Moreover, 
this maximum is globaoif T{E) decreases with energy, see Eq. ( I23i) . Thus, at pc,i < Po < Pc,2 
the integral for the transmission probability Voo is saturated at Ep = Vq, i.e. precisely in the 
vicinity of unstable semiclassical trajectory. 

The above considerations deserve two remarks. First, one can show that the semiclassical 
trajectory saturating the probability ( |22i) satisfies the boundary conditions of the previous 
Section, Eqs. ([9]), ( ITOl) . This is natural, since Eqs. ([9]), ( JTOl) are obtained from the requirement 
that the probability is extremal, see Appendix |Al In Appendix [C] we explicitly solve the 
semiclassical equations in the model (|2T|) and find the trajectory x{t). It is unstable if 

Pc,l <P0< Pc,2- 

^The case of "pathological" potentials with non-monotonic T{E) should be considered separately. 
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Second, the interval of unstable motions is semiclassically large, pc,2 — Pc,i ~ ^/g, and 
shrinks to a point as 0"^ — 0. This means that unstable semiclassical dynamics is essential 
in the case (7^ <C 1, cr^ ~ 1 which we consider. 

Let us compute p(r) in the case of unstable semiclassical trajectories x(t). We have 
learned that such trajectories do not interpolate between the in- and out- regions of the 
process, but start in the in-region at t = and approach exponentially the barrier top x = 
as t —>■ +00. At large times one writes x{t) = c_e~^"*, where c_ is a complex constant. 
Following the method of the previous Section, we introduce modification, Eq. f|T5l) . We will 
restrict attention to the case of small e which will be important for calculating p(r). One 
notes that the modified trajectory x^{t) moves in the comp/ex potential 

U,{x)=U{x)-^-e{x2-x) . (25) 



It cannot stop at the barrier top due to conservation of energy, which is real by Eqs. (fTOj) . 
The evolution described by x^(t) proceeds in three stages: the trajectory arrives into the 
region x ~ 0, spends some time there and then departs for the out-region. The modification 
term can be neglected during the first and last stages of evolution, but destroys unstable 
motion near the barrier top. Indeed, at Xe ~ one writes, 

x,{t) = c_ e-^-* + ec+ e+'^'* , (26) 

where e in the second term reflects the fact that the original trajectory is recovered in the 
limit e — ^ 0. Since e <C 1, the time of motion near the barrier top is large, and one computes 
Tint, Eq- ( |T6l) . using the evolution ( l26l) . 

r = -— ln(e/eo) + 0(6), (27) 

UJ- 

where eo ~ 0{g^) is some constant. In practice one finds the value of this constant explicitly, 
by obtaining the trajectory at all three stages of the modified evolution and calculating 
Tint[X^. Given e(r), one computes the suppression exponent from Eq. f[T5]) . 

r" + 00 



r 2e 
= + 2 / cir'e(r') = F^ + ^i 

Jt ^- 



where F^o is independent of r. 

Now, we estimate the modified prefactor aX e <^ 1. The procedure for evaluating 
is described in Appendix [Cl we proceed by following this procedure in the case of small e. 
One finds the perturbation which satisfies the linearized equations of motion in the 

background of the modified trajectory Xeit). Boundary conditions for ipeif) are imposed 
in the out-region. Evolving the perturbation back in time, one notes that V'e(^) grows 
exponentially during the second stage of the modified evolution, Eq. (1261) . Thus, it becomes 
exponentially large at t = 0, ijjei^) ~ exp(co'_r). At e -C 1 the prefactor is estimatecl^ as 
A, ~ [?/',(0)]-^/^ see Eq. ([63]) of AppendixO One finds A, ~ exp(-cj_r/2). 



^°We disregarded irrelevant numerical factors in Eq. ([55]) and perturbation Xe(t) which does not grow 
exponentially. 
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Substituting e(r), and into Eq. (fT7|) . we obtain the distribution 

2eo 

- exp \ —LU^T 



Pu{r) 



(28) 



which has the form of Gumbel distribution of type I [35]. We will argue below that Eq. (I28l) 
is valid in the region r ^ | logeo|/cc;_ where Pu{t) is not exponentially small. By using the 
subscript u in the notation p«(r) we stress that Eq. fl28|) is valid in the case of unstable 
semiclassical trajectories. 

Several remarks are in order. First, Pm(t) is drastically different from the respective 
Gaussian distribution flTIJl) which is valid in the case of stable semiclassical trajectories. 
Indeed, function (l28l) is essentially asymmetric. Besides, it corresponds to large mean time 
of transmission and large time dispersion. 



-In 



260 



+ 



(29) 



where 7 is the Euler constant. The quantities (!29l) depend on the semiclassical parameter 
as (r) ~ I log^f^l, cr^ ~ 0((y'°), in contrast to the respective scalings (r) ~ 0{g^), a"^ ~ g"^ in 
the case of stable trajectories, cf. Eqs. (|2U]) . 

Second, the form of the distribution (125]) is universal. Indeed, in deriving p„(r) we did 
not use any information about the potential U{x) or initial wave packet '^i{x). All this 
information is encoded in two parameters of the distribution, eo and uj-. Moreover, one 
represents Eq. (1281) in the form 



Pu{t) = id- exp {-t^-(r - (r)) - 7 



-<^_(r-(T))-7 



} . 



where the first of Eqs. fl29|) was used. All parameters in this formula disappear after appro- 
priate shifts and rescalings of r. 

Third, during the derivation of Eq. fl28l) we assumed that e ^ 1. This is legitimate at 
r ^ |logeo|/a;_, see Eq. (|27|) . In particular, the region Pu{t) ~ 1 corresponds to r ~ (r) 
which leads to e ~ g"^, cf. Eqs. ( 129]) and (127]) . 

We derived expression (125]) using the properties of the modified semiclassical dynamics. 
This expression is illustrated in Appendix [C] where we compute p(r) for the mode l (|2T]) . In 
the region pc,i < Po < Pc,2 the result of Appendix [C] coincides with Eq. (125]) . wherqlj 



eo 



sm 



gvQ 



(30) 



Pc,i and p = Pc,2 = V^/g, see Eq. I^M)- Function (125]) is depicted 



Note that eo vanishes at p 
in Fig. Hi, sohd hnS 

^^Exponential dependence of eg on X2 — xi has simple physical meaning. Namely, Eqs. (j29p . (|30p imply 
that (r) = g(a;2 — xi)/ \/2 + . . . , in accordance with the fact that the modes saturating the total transmission 
probability move with momentum p « V2Vo = V^/g. 

^^We plot the distributions as functions of r — (r); in this case all parameters except for a;_ vanish from 
Eq. (|^5]). Besides, we use lu^ = \/2 which is valid for the potential (|2ip . 
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Figure 3: (Color online) (a) One-dimensional traversal-time distributions in the case of 
unstable semiclassical dynamics. Points and solid line represent exact quantum mechanical 
results and Eq. fl28l) . respectively, (b) Average time of passing in the model i^; = 0.002, 
a2 = 0.02, gx2 = -gx, = 10. 

As a check of the modified semiclassical method, we compare Pn(T) with exact quan- 
tum mechanical results (points in Fig. [3^). The latter results are obtained by propagating 
numerically!^ wave packets in the model fl2T|) and using Eqs. ([1]), ([2]). One sees that the 
exact graphs are almost symmetric at large g^ but change their form and approach Pm(t) as 



Our semiclassical expressions for the average time of passing are summarized in Fig. [3b 
which displays dependence of (r) on pq in the model ( 12T|) . The interval of unstable semi- 
classical motions Pc,i < Po < Pc,2 is delimited on the graph by the vertical dotted lines. The 
regions to the left and to the right of this interval correspond to stable tunneling trajectories 
and classical over-barrier solutions, respectively. In accordance with the above discussion, 
we compute (r) by Eq. (12^ in the central part of the plot and by Eq. (12(1 in the rest of it. 
The corresponding dependencies are marked in Fig. [3)d as "unstable" , "stable" and "classi- 
cal". All three graphs exhibit unphysical behavior in the vicinities of critical momenta pc,i, 
Pc,2, where the "stable" and "classical" curves grow to infinity, and "unstable" graph sharply 
drops down. This behavior is related to "phase transition" which transforms stable trajec- 
tories into unstable and Gaussian distribution ([TI?]) into Eq. (1251) . Near the critical points 
(r) is computed with the full semiclassical distribution ([T7|) . This calculation is performed 
in Appendix O In Fig. [3]d we show the result by the solid line which smoothly interpolates 
between the "stable" , "unstable" and "classical" graphs and coincides with the average time 
of passing extracted from the exact quantum mechanical calculations (points). 

-"^^This can be achieved by Fourier methods [36] or by using the basis of energy eigenstates in the model pip . 
We checked that both methods produce equivalent results. 
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5.2 Two dimensions 



We proceed by considering the case of sphaleron-driven transitions in the two-dimensional 
model of Sec. [3l At the level of complex trajectories these transitions look similar to one- 
dimensional activation processes of the previous Section. However, there is one important 
difference: the origin of semiclassical instabilities in the sphaleron-driven case is related 
to non-linear interaction between the degrees of freedom, rather than to the momentum 
dispersion in the initial state. Recall that the new mechanism is relevant in the region 
Pc,i < p < Pc,2, where Pc,i and Pc,2 are critical momenta defined in Sec. [31 

In the semiclassical calculations of this Section we use numerical methods of Refs. 
Namely, we smoothen the step-function in the modified potential, Eq. 

9{x) ^ 9,{x) = , 

where the width a is supposed to be smal0. After that we introduce non-uniform lattice 
t G {tfcl k = 1 . . . Nk}. The trajectory Xg(t) is found from the modified equations of motion 
and boundary conditions IQ, f fTOl) in the discrete case. Then, the values of F^, and p(r) 
are computed by Eqs. ( 137]) . ( l42l) of Appendix |A] and Eq. ( fTTl) . respectively. 

We start with the simplified case when the total energy E of the original sphaleron-driven 
trajectory x()f:) is close to the height of the barrier, — Vq -C Vq. The sphalerons with these 
energies describe small linear oscillations around the saddle point x = y = of the potential; 
for the modified semiclassical motion in the vicinity of the saddle point one writes, 

x,(t) = x_ (c.e""^-* + ec+e""-*) + x+ccos(cj+t + ^) , (31) 

cf. Eq. (1261) . Vectors x_|_ and x_ in this expression run along "stable" and "unstable" 
directions of the saddle point, while uj± represent respective "frequencies." Using the same 
arguments as in the previous Section, one shows that the distribution p(r) is given by Eq. (p8l) 
if Eq. (EH) is justified (i.e. at E - K) < Vq). 

In Fig. m we compare one-dimensional formulc0, Eq. ( |28ll . with two-dimensional semi- 
classical distributions ffTTj) at different values of g^. Note that all distributions in Fig. [Hare 
semiclassical; nevertheless, their forms depend on g"^. One explains this as follows. We se- 
lected parameters of Fig. H] in such a way that the energy E of the original trajectory (e = 0) 
is almost equal to the barrier height, {E — Vq)/Vo ~ 10^^. On the other hand, the energies 
of X(:(t) remain close to Vq at small e (large r) and depart from it at e ~ 1 (smaller r). Thus, 
linear approximation (l3T|l breaks down if the value of r is not sufficiently large. Due to this 
fact the semiclassical distributions in Fig. H] coincide with Pu{t) only at small g^. Indeed the 
central parts of the distributions correspond to r ~ (r) ~ | log^f^l, see Eq. ( 129|1 : at large g"^ 
the distributions shifS to smaller r, where Eq. (l3Tll gets violated. 

-"^■^Practical calculations show that p(t) is almost independent of a: at a = 0.1 the values of p stabilize at 
the level of accuracy 10~^. In what follows we use a = 0.1. 

^^We find Pu{t) using the "negative" eigencurvature = 1 — + ^/l + uA/A of the potential (O at 
the saddle point. 

^^This is not seen in Fig. 0] where the distributions are functions of r — (r). 
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Figure 4: (Color online) Comparison between one-dimensional semiclassical formula (!28l) 
(solid line) and semiclassical traversal-time distributions in the case of splialeron-driven 
tunneling (dashed lines). The parameters of the graphs are gpo = 1.038, g^Ey = 0.45; 
= 0.005, gx2 = -gxi = 20. 

In general case when the energy E of the original unstable trajectory x(t) is substantially 
higher than Vq the distribution p(r) is not given by Eq. ( l28l) . regardless of whether g^ is small 
or not. However, many qualitative features of Pm(t) are valid in the higher-energy case as 
well. Consider e.g. Fig. [5^ which corresponds to {E — Vo)/Vo ~ 0.2. The distributions 
in this figure are highly asymmetric; they have long exponential tails at large r and steep 
front-ends. Exponential behavior of p(r) at large r is illustrated in Fig. [5]d where Fig. [5^ is 
replotted in logarithmic scale. One writes, 

p(r) oc exp(— const ■ r) , (32) 

since the graphs in Fig. [5}d are almost linear at r ^ (r). 

In Figs. [5] we compare the semiclassical (lines) and exact quantum-mechanical (points) 
results for the traversal-time distributions; one observes agreement. The exact results are 
obtained by evolving numerically wave packets in full quantum system and implementing 
Eqs. ([1]), ([2]); we describe the respective numerical method in Appendix [Dl 

The above properties of p(r) naturally lead to large mean times of transmission (r) 
and large time dispersions cr^ in the case of sphaleron-driven tunneling. In Figs. [6] we 
plot (r) and cr^ as functions of g"^ at different values of Pq. One observes distinct behavior 
in the cases po < Pc,i and po > Pc,i corresponding to stable and unstable semiclassical 
dynamics, respectively. At small po the graphs in Figs. [6] agree with Eqs. fl20l) : (r) remains 
constant, and cr^ vanishes in the limit g^ —>■ 0. In the unstable case po > pc,i one sees that 
(r) ~ I log I while cr^ grows as g^ — > 0. Using these dependencies, one can discriminate 
between the standard and sphaleron-driven mechanisms of tunneling. Note that analogy 
with one-dimensional activation processes suggests logarithmic growth of (r) but constant 
cr^, see Eqs. ( |29|) . Apparently, the unexpected behavior of cr^ in Fig. [6)d is related to features 
of non-linear multidimensional dynamics which deserve separate study. 
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Figure 5: (Color online) (a) Semiclassical (lines) and exact quantum-mechanical (points) 
traversal-time distributions in the model ([Tj). The distributions are plotted in units of their 
maximal values Pmax- Numbers near the graphs give values of g"^; the other parameters 
are gpo = 1.52, g'^Ey = 0.05; (7^ = 0.005, gx2 = —gxi = 20. (b) The same as (a) but in 
logarithmic scale. 
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Figure 6: (Color online) Mean times of passing (a) and time dispersions (b) as functions of 
g'^ in the model ([7]). Numbers near the graphs give values of gpo- The other parameters are 
the same as in Figs. [5l Note that gpc,i ~ 1.421; graphs with po > pc,i correspond to the case 
of sphaleron-driven tunneling. 
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Figure 7: (Color online) Mean times of passing (r) and rescaled time dispersions o^jg^ as 
functions of g-p^. The semiclassical and exact quantum-mechanical results are represented 
by lines and points, respectively. The parameters of the graphs are the same as in Figs. O 
The critical value g'Pc,\ ~ 1.421 is shown by the vertical dotted line. 

In Figs. [7] we plot (r) and o"^! g^ as functions of mean initial momentum po- Dashed lines 
in these figures represent quantities computed with the distribution (fT7|) at different values 
of g^-^ they agree with the respective exact quantum-mechanical data (points). One sees that 
all graphs in Figs. [7] coincide in the region po < Pc,i corresponding to stable trajectories. This 
is what one expects, since in stable case (r) and d^jg^ are given by Eqs. (!20l) and therefore 
do not depend on g^ . "Stable" expression^ (l20l) are represented by the solid lines. Above 
the critical momentum pc,i the situation in Figs. [7] changes drastically: both (r) and o'^jg^ 
depend on the semiclassical parameter, so that small g^ correspond to large mean times of 
passing and large time dispersions. 

Finally, let us repeat that the processes of sphaleron-driven tunneling survive in the 
limit o"p — > 0, in contrast to one-dimensional activation processes of the previous Section. 
Our numerical calculations show that, indeed, at arbitrarily small there exists a region 
of unstable semiclassical motions pc,i < Po < Pc,2; the width of this region stays finite in the 
hmit cTp 0. 

5.3 Long— time behavior of traversal— time distribution 

Let us return to the one-dimensional setup of Sec. [3] and consider transmission at high initial 
momentum > V2Vo. This process is described by stable over-barrier trajectory xii) which 
passes the region xi < a; < 0:2 in finite time r = (r). In the vicinity of (r) the distribution 
p(r) is Gaussian, see Sec. HI However, we are interested in the region r — * +00 where the 
value of p(r) is exponentially small. We are going to show that in this region p(r) does not 

""^ ''Apparent singularities of "stable" graphs at po — * Vc.i indicate that Eqs. (|20|) break down in the imme- 
diate vicinity of pc,i- 
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Figure 8: (Color online) Distribution p(r) plotted for over-barrier transitions in the 
model ([21]); = 0.002, po = lAh/g > ^JWq, = 0.02, gxi = -gx2 = 10. 

resemble the Gaussian hat at all; rather, its form is similar to Eq. ( |28l) . 

Semiclassically, the asymptotics of p(r) is elucidated by studying the properties of the 
modified trajectories x^{t). At large r these trajectories linger near the barrier top, see 
Eq. fl26|) . One uses considerations of Sec. 15.11 and arrives at the formula 



cf. Eq. ( I28l) . This asymptotics is universal, since parameters U-, eo and poo disappear after 
appropriate shifts and rescalings of r and p. 

Let us point out one important difference between Eqs. (I33l) and fl28|) . In Eq. fl33l) we 
consider the case of over-barrier transmissions when p(r) reaches its maximum at finite 
T = (r). As one goes away from the maximum, the exponent in Eq. f|T7|) grows. Since 
r-derivative of this exponent is equal to — 2e, the parameter e is negative at r > (r); thus, 
eo < by Eq. fl27|) . One concludes that the asymptotics fl5^ monotonically decreases with 
T and Poo is exponentially small. 

In Fig. [8] we compare Eq. (!33|) (dashed line) with full semiclassical distribution (ITTl) (solid 
line) and exact quantum mechanical results (points). One observes coincidence at r (r). 

6 Summary and remarks 

The qualitative result of this paper is summarized as follows: quantum transitions take 
long times if they are described by unstable semiclassical trajectories. To understand how 
long this "long" is, we introduced the semiclassical parameter g"^ in Eq. (jl]). Then, unstable 




(33) 
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semiclassical dynamics can be identified in the limit ^ where one obtains logarithmically 
large mean times of transitions and large time dispersions. This behavior is drastically 
different from that in the case of stable trajectories. 

We defined traversal time from time-of-arrival measurements by the asymptotically dis- 
tant detector. One may be unsatisfied with this definition and use another one, see e.g. 
Refs. [m |12l [131 [IS [IS] ■ However, the qualitative conclusion about large traversal times 
should be valid for any reasonable traversal-time definition, since large time delay comes 
from the almost-classical motion in the vicinity of the unstable mediator orbit (sphaleron 
or the barrier top). 

The starting point of our analysis was general semiclassical expression f[T71) for the 
traversal-time distribution p(t). In Sec. [S] we applied this expression to three processes 
related to unstable semiclassical dynamics: one-dimensional activation transitions, multidi- 
mensional processes of sphaleron-driven tunneling and over-barrier transmissions with large 
delays. The results were obtained analytically in one dimension, Eqs. (128!) and (!33|) . and 
numerically in the case of sphaleron-driven tunneling. We observed the following feature due 
to semiclassical instabilities: p(r) is highly asymmetric and contains long exponential tail 
stretched toward r — >• +oc. This distribution is drastically different from Eq. f[T9|) obtained 
in the case of stable trajectories. 

Note that our results are somewhat different from those of Ref. [9] where power-law de- 
pendence of |\E'(r, a;2)P on r was reported. This dependence should be compared with the 
exponential tail of p(r) in Fig. \5jp. It seems that the difference can be explained in the fol- 
lowing way. In the model of Ref. |9] many semiclassical trajectories contribute substantially 
into \l/(r, X2); one assumes that the number of these contributions grows with time r. On the 
other hand, only one trajectory is relevant in our model (jTj). Summation over large number 
of semiclassical contributions may be responsible for the modification of large-r behavior of 
wave function in the model of Ref. [5] as compared to our case. 

Unusual form of p(r) and related scalings of (r), with can be used for experimental 
identification of sphaleron-driven tunneling and one-dimensional activation transitions. 
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A Semiclassical transmission probability 

Let us evaluate semiclassically the total probability Voo- For definiteness we consider two- 
dimensional model of Sec. [3l semiclassical calculations in one dimension are performed in 
the same way (see comments in Appendix [Cll. 
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One starts with the path integral for the final state, Eq. ([8]). The in-state entering 
this integral is described in Sec. [31 In the semiclassical approximation one writes, 

vi>.(x,) = . e^^(-"^»)+-/4 , (34) 

(2vr)3/^V^ 



where Py{yi) = ^/'^Ey — oo'^y'^ is ^/-component of the initial momentum, and 

/Vi 
Py{y')dy' , (35) 

is in-state correction to the classical action. 

Note that the integrand in Eq. ([8]) contains fast-oscillating exponent exp[z(S' + -B)], 
where B comes from Following the prescription of the saddle-point method, one finds 
the trajectory x(t) extremizing the functional S + B. In particular, x(t) should satisfy the 
classical equations of motion and Eqs. 0. After the saddle-point integration one obtains, 

^(r, x^) = D-^'^ ■ e^(^+^)+-/4 ^ (36) 

where all prefactors and saddle-point determinants are collected in Di\ we evaluate them 
below. We substitute Eq. fl36l) into the probability formula ([1]) and compute the saddle- 
point integral over xj. In this way we derive the conditions ffTOl) and probability expression 
f fTTj) . [Note that r is sent to infinity in the final formulas.] The suppression exponent of the 
probability is given by the value of the functional 

F = 2/lm(5 + 5) (37) 

on the trajectory x(t). The probability prefactor 

A = r-^^ (38) 



\DIWD2 

involves additional determinant D2 due to the saddle-point integration in Eq. ([1]). 

Before proceeding with the factors Di and D2, we preview the result for the probability 
prefactor A. One considers the set of perturbations 5x(t) satisfying the linearized equations 
of motion 

5i + f/"(x(t))5x = (39) 

in the background of complex trajectory x(t). The basis in this set is formed by the "mo- 
mentum" and "coordinate" perturbations il)^^\t) and ^^-'■'(t), j = 1,2. We assume that the 
basal perturbations satisfy two conditions: (i) they are normalized by 

where f2 = (ip A rfx is the canonical symplectic form; (ii) they are real in the out-region. 
Note that both of the above conditions can be imposed in the out-region and therefore 
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constitute final Cauchy data for tlie perturbations il)^^\t) and X^^\'^)- Note also that the 
basal perturbations are not completely specified by (i), (ii); in practice one fixes the freedom 
by supplying some additional data. 
One introduces the functionals 



(5/i[(5x] = 5xi 



(40) 



which measure changes of initial conditions dHj) caused by the perturbations (5x(t). Changes 
due to the basal perturbations form two matrices, 



The probability prefactor is written in terms of these matrices, 

A = 



^27r-det[zen+-zne- 



(41) 



(42) 



This expression is valid for any set of basal perturbations i/j^-^-*, x^-') satisfying (i), (ii). 

Semiclassical calculation of transmission probability Voo is summarized as follows. One 
finds the complex trajectory x(t) from the classical equations of motion and boundary condi- 
tions ([9]), ffTOl) . The exponent F is given by the value of the functional fl371) on this trajectory. 
One also constructs four linear perturbations "y^^^ by solving Eq. ( l39i) with the Cauchy 
data (i), (ii) in the out-region. Prefactor A is computed by Eq. ( l42l) : it involves changes 
of initial conditions due to the perturbations i/?^-''' and x^-'-'. Given A and F, one finds the 
probability Poo from Eq. (ITTl) . 

We derive Eq. fH2|) in two steps. First, we completely specify the basal perturbations 
ijj^^\t), X^'^K't) ciiid find A in terms of these perturbations. Second, we prove that Eq. 
is basis-independent once the conditions (i), (ii) are met. 

One computes Di by collecting the prefactors in Eq. ([8]), 



Di = {27T 



1 3/2 Vi 

2uja,r 



X {2ni) 



det 



9xy 

9x,; 



{27Ti) 



det 



d\S + B) 
9x? 



(2vr) 



3/2 Vi 



■det 



d\S + B) 



2uoap dxfdxi 



(43) 



where xj is kept constant in the differentiations. Three factors in the first line come from 
the initial state and saddle-point integrals ovei0 x(t) and Xj, respectively. Evaluating 
explicitly the derivatives of S and B, one obtains. 



{2n 



1 3/2 



2ujar, 



det 



1 

m 



dii-f 



-2ial 



\ d^i 



(44) 



=const 



^We performed integration over x(i) using the semiclassical Van Vleck formula [57] . 
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where we wrote determinant of diagonal matrix and performed matrix multiplica- 

tion. We choose the basal perturbations i/y^^\t) = dyi(t)/dxf, ijj^'^\t) = d:s.(t)/dyf, where 
the derivatives are taken at fixed xj. Clearly, these perturbations satisfy the linearized 
equations (1391) . Their final Cauchy data are 



V,«(r) = i/'^'^(r) = 0, 
In terms of i/^-perturbations 



(r) 



(27r) 



3/2 



detn 



(r) 



(45) 



where 11 is 2 x 2 matrix defined in Eq. fHT]) . 

Determinant D2 comes from the integration in Eq. ([T]), 



Do = det 



.d±f 

-t-:r-^ + h.C 



(9x 



/ 



(46) 



where the derivatives are taken in the subclass of trajectories satisfying Eqs. (jnD- Let us define 
two remaining basal perturbations as X^^K't) = dx(t)/dxf, X^'^K'^) — dxi^)/dyf, where xj is 
kept fixed in the differentiations. These perturbations satisfy 



X^'\r) = x'Hr 



(2), 



. 



(47) 



Note that the matrix in D2 is not directly related to ^-perturbations, since the differen- 
tiations in Eq. fH6l) are performed with Eqs. ([9]) kept fixed. One introduces two auxiliary 
perturbations 



/^«(t) 



9x(t) 



Eqs. |9l 



dyf 



(4J 



Eqs. |9j 



which enter Eq. fl46|) . They can be decomposed in the if), x-basis, 

2 



(49) 



k=l 



where we took into account Eqs. (H5|) . (H7!) . Since K-perturbations do not change the initial 
conditions, one writes SIjIk^''^] = 0. Substituting representation (11^ into these equations, 
one obtains 

= -n-'Q . 

We finally write D2 in terms of the perturbations k,^^^ and use Eq. fH^ . We find, 

D2 = I det n|-2 . det[ien+ - in0+] . 



26 



Using expressions for Di, D2 and Eq. (!38|) . one gets Eq. (Il2i) . 

We finish this Appendix by proving that any set of basal perturbations with properties 

(i) , (ii) can be used in Eq. ([42]). Note first that V^^^ X^^'' introduced in Eqs. do 

satisfy these properties. Denote by x*-*^ some other set which also passes (i), (ii). One 
decomposes new perturbations in the old basis, substitutes them into Eq. (H2!) and uses (i), 

(ii) . In this way one proves that the value of the prefactor (H2!) does not depend on the 
choice for basal perturbations once they are real in the out-region and correctly normalized. 



B Expression for p{t 

Here we give details of the semiclassical evaluation of p(r) which was outlined in Sec. |H We 
start by computing the integral in brackets in Eq. (fT4l) . The result is (cf. Eq. (!36l) ) 



^{t, X.) = r dr' ^ D-y^ ■ e^(5.[x.]+B[x.])+iV4 ^ (50) 

Jo Jioo '^T^9 

where 8^ = 3 + ie(Tint — T')/g'^ is the modified action. By subscripts e here we mark the 
quantities computed with the modified trajectory Xe(t) and modified action S^. Recall that 
Xe(t) extremizes 5*^, arrives at Xe(r) = and satisfies the initial conditions ([9]). 

Now, we evaluate the integral over e in Eq. (!50|) . Since Xe(t) extremizes + B, one 
writes, 

^{S, + B) = ^- ^{S, + B) + |-(S, + B)= ^(T,„, - r')/g' . (51) 
ae ae ox^ oe 

Thus, the exponent in Eq. fl50l) is extremal with respect to e when Tj„t[xe] = r'. Performing 
the saddle-point integration over e, one obtains. 



v&(r, X,) = r -^=J~^ ■ e^(^^+^)+-/^ . (52) 
Jo g^2^d\J dr' 

Additional prefactor in this expression arises due to fixation of Tint- 

Next, we substitute \E'/ into the probability formula ([1]). Note that Vt involves, besides 
\E'(r, xj), its complex conjugate. Starting with the path integral for the latter, one shows 
that expression for ^*(r, xj) is obtained from Eq. flS^ by changing the signs of the exponent 
and e. We write, 



Vr= I d^f r ^ r-^^ /jl^ . ^.(5.[x.]+B[x.])-.(5_.[x_.]+B[x_.]) (53) 

where f and e come from \E'*(r, xj); they are related by Tj„j[x_f] = r. 

One sees that the probability formula involves integrals over two traversal times r' and 
f spent by the respective trajectories in the region x < X2- We introduce new variables 
r_|_ = (r' + f)/2 and r_ = t' — f and perform the integral over r_ in the saddle-point 
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approximation. This leads to the condition e = e. Note that after the integration the value 
of e is related to r+ by the implicit relation 

T,„4[x,]+T,„4[x_,] =2r+ . (54) 

Assuming that the trajectory Xe(t) is unique, one proves that the saddle-point value of e is 
real. Indeed, it is straightforward to check that x* = x_e*- Thus, for real e the function 
r_|_(e) is real, see Eq. 0541) . The same is correct for the inverse function e(r_(_). Equation 0541) 
and reality of e lead to the relation 

ReTj„4x,] = r+ (55) 

which implies that r_|_ is the real part of the time interval spent by the modified trajectory 
in the region between xi and X2- 

Finally, let us evaluate the integral over xj in Eq. fl53|) . It is taken in the same way as in 
Appendix |Xl However, we obtain conditions at finite r, 

ImXe(r) = ImXe(r) = (56) 

cf. Eqs. (do]). 

We arrive at the following representation. 



de 



dr. 



A, ■ e"^^/^' , (57) 



where only the integral over r+ is left. Quantities and in this formula are given by the 
same expressions (137|) and (H2|) . but with the modified action (fT5|) and modified trajectory 
Xe{t). Note that now r+ enters the modified action instead of r'. 

In Refs. [2S1 representation similar to Eq. (IFr|) was used for evaluation of total transmis- 
sion probability Voo- In this case one sends r to infinity and performs the integral over r+. 
We follow another path and calculate p(r) by differentiating Eq. fl57|) with respect to r. 

Let us show that the integrand in Eq. fl57|) does not depend on the duration of the 
process r. Consider first the modified trajectory Xe(t). Parameter r enters the boundary 
value problem for Xe(t) via the conditions fISBl) imposed at t = r. Note, however, that 
Xe(r) lies inside the detection region x > X2. Indeed, due to Eq. r+ is equal to the 
time of motion in the region x < X2; at r+ < r the trajectory leaves this region before 
t = T. From Eqs. (fT5l) . (fT2l) one finds that the classical equations of motion are not modified 
in the detection region x > X2. Thus, Xe(t) describes real classical evolution at t ~ r, 
and conditions of reality, Eqs. (156!) . can be imposed at any point of this evolution without 
affecting the trajectory. One concludes that Xe(t) does not depend on r. The suppression 
exponent F^, Eq. 0371) . is independent of r as well, since it involves the imaginary part of 
Se computed on the modified trajectory. Function e(r+) is defined in Eq. f l55|) : it does not 
depend on r due to definition of Tint, see Eq. f|T2l) . 
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Now, we change the duration of the process, r — >■ r + 6t, and find the respective changes 
of the basal perturbations X- For concreteness we use the perturbations defined by 
Eqs. ( H5|l . (H71) . Solving explicitly Eqs. (!39|) in the asymptotic region x > X2, we find that 

A -5r\ /t/^^'^A f cos{uj6t) -sin(cu5r)M /V^'^A 

U^'V " VO 1 J W'U ' U^'V ^ Wsin(^5r) cos(^5r) J [x^^) J ' 

(58) 

due to 6t. Using these transformations, one exphcitly proves that A^, Eq. fH2|) . does not 
depend on r. 

Now, take a look at the probability Vr, Eq. (!57|) . We have shown that the quantities F^, 
e(r+), in the integrand of this expression do not depend on r. Thus, the derivative p(r), 
Eq. ([2]), is simply given by the integrand at r_|_ = r. In this way one gets the semiclassical 
expression (fT7|) of Sec. IH the relation (ITB]) is obtained from Eq. (1531) . 

Finally, we prove Eq. ffTSj) : 

since is extremal with respect to Xe(t) and e. 



C Semiclassical calculations in one dimension 



In this Appendix we perform explicit semiclassical calculations in the model (I2T1) . 

For a start, we find the original complex trajectory x{t). General classical solution with 
energy E < Vq = l/g"^ looks like 



sh{gx) = ch gV2E{t - t 



(60) 



where E and to are real by Eqs. (fTOj) . One notes that the solution (160!) . if taken along 
the real time axis, describes refiection from the potential barrier; trajectory with correct 
asymptotics is obtained by winding the time contour around the nearest upper branching 
point of the solution. From the initial condition ([9]) one findS the energy of the solution, 
E = {po + 27rcrp/ (7)^/2. Since E does not exceed the height of the potential barrier, the value 
of po is smaller than p^^i, Eq. flM|) . 

The case E > Vq corresponds to classical over-barrier motions which trivially produce 
E = pI/2 and therefore Pq > \/2Vo = pc,2- One concludes that transitions in the region 
Pc,i < Po < Pc,2 are not described by the trajectories with E 7^ Vq. The remaining solution 



sh.{gx) = — exp —V2{t — to) 



(61) 



^^Euclidean parts of this contour correspond to motions under the potential barrier. 

^°Recall that initial conditions are imposed in the asymptotic region x —00 where motion is linear. 
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corresponds to E = Vq = l/g"^. It starts in the in-region and approaches the barrier top 
as t ^ +00. Note that Eq. (!6T!) automatically passes the final boundary conditions (fTOj) . 
Thus, to is complex; one finds it from Eq. Qj. Trajectories (16T!) exist at any pq. Calculating 
the suppression exponent (jSZD, one explicitly checks that they are sub-dominant outside the 
interval Pc,i < Po < Pc,2 and describe transitions through the barrier if pq belongs to this 
interval. 

Note that the fixed-energy transmission coefficients Vp are exactly known in the model 
f l2T|) . Using these coefficients and Eq. fl22|) . one explicitly confirms the semiclassical expres- 
sions for Voo, Refs. [26], in cases of stable and unstable semiclassical dynamics. 

We proceed by introducing modification. One notices that the modified potential fl2Bl) 
involves step-function and therefore has two analytic continuations, which start from the 
regions x < X2 and x > X2. We find the corresponding parts of the modified trajectory and 
sew them and their time derivatives at t = ^2 where Xe{t2) = X2- Note that the energy E of 
trajectory is conserved across the sewing point. The time ^2 and energy E are real due to 
the final boundary conditions fllOp . Moreover, ^2 = ^ by Eq. f[T^ . At t < ^2 the modification 
adds the imaginary constant —ie/g"^ to the potential. Thus, in this region x^{t) is given 
by Eq. fl60|) where one should substitute E ^ E^ = E + ie/g"^. Note that parameter to is 
complex after modification: it is found from Xe{T) = X2. Using this condition and Eq. ([9]), 
one obtains non-linear complex equation 

f{Ee) = ^ [V2E, - Po) +X2-XI- + ^ • In ^^f^y = (62) 

which relates E and e to po and r. 

Now, we compute p(r) using the modified quantities e, F^, and Eq. (|T71) . The function 
e(r) is found from Eq. (|62l) . The modified exponent F^ is computed by substituting Se and 
x^(t) into Eq. fl37j) and performing integration. [We disregard the y-dependent part of the 
functional B, Eq. fl35l) .] The final ingredient is the probability prefactor A^ which is given 
by one-dimensional formula 

A = (63) 

^imi5inx]-5hm ■ 

Here by omitting the subscript e we mean that Eq. (163!) can be used both in modified and 
unmodified cases. Note that one-dimensional prefactor can be obtained from Eq. (l42l) by 
reducing the size of matrices 11 and G and omitting the kinematical factor cu / ^/2tt. 

Linear perturbations Xe, "^e are calculated by differentiating the trajectory Xe{t) with 
respect to parameters. 



/2E ' ^ ' dpo ' 

where we have already used the Cauchy data (147|) for x(t). Conditions (145!) give (3~^ = 
dV2E/dpo, where the derivative is taken at e = const. Substituting the perturbations into 
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Eq. (!63|) . and using Eq. ([9]), one arrives at the expression 

for the probabihty prefactor. 

We substitute and into Eq. (fTTl) and use some algebra. The result is 

where Af' is normalization factor. Note that f{E^) is defined in Eq. ( 1621) . and the values of 
i?, e are found from /(-EJ = 0. 

The expectation value (r) of the distribution fl64|) is plotted in Fig.[3)D, solid line. "Stable" 
and "unstable" distributions, Eqs. (fT9|) and (128|1 . (l30l) are obtained by considering Eqs. (162|) . 
in the region of small e. 



D Exact results in two dimensions 

We compute exact traversal-time distributions in the model ([7]) using the following numerical 
procedure. 

First of all, we introduce the basis of stationary eigenfunctions 'il)p{x,y) satisfying time- 
independent Schrodinger equation with scattering boundary conditions, 

y) — i> ipyiii) ■ e'^^ + reflected waves , as x — — oo , , . 

y) —i> outgoing waves , as x ^ +oo , 

where '^/'^(y) is the initial oscillator state with energy Ey. The total energy corresponding to 
iljp{x,y) is Ep = p'^/2 + Ey. 

One decomposes the initial wave packet (Sec. [2D in the above basis. 



*i(a;, y) = i/j^{x) ■ ^/jy{y) = -j= / dp i/j^{p) e'^"" ■ ipyiy) = — = / dp ip^{p) ^/jp{x, y) , 

\J I'K J y 111 J 

where we performed Fourier transformation and used the asymptotics (!65l) of ipp{x, y). Then, 
we propagate the wave packet from t = to t = r. 



*(r, X, y) = —= / dp i)^{p) i)p{x, y) ■ e '"^^ . (66) 
y 2,71 J 

Given representation 0661) one computes the total probability current through the line x = X2- 
The distribution p(r) is proportional to this current, see Sec.[2l the proportionality coefficient 
is obtained from the normalization condition J p(r) dr = 1. 
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In practice we computed the exact distributions p(r) alF^ = 0.1 and 0.01. At each value 
of we generated the set of eigenfunctions ippi^yV) by solving numerically the stationary 
Schrodinger equation at several energiej^ Ep. To this end we used the numerical method of 
Ref. [33] (see Ref. [M] for the Fortran 90 code). Given the eigenfunctions, we constructed 
the linear combination fl66p and calculated the total probability current through the line 
X = X2- Computing numerically the normalization integral, we finally obtained p(t). 

There are two limitations of our numerical procedure; both of them are relevant for the 
numerical data at g"^ = 0.01. First, our method of solving the Schrodinger equation fails 
to produce correct ei genf unct ion at g"^ = 0.01 and Ep < 1. Accordingly, the exact data 
at g"^ = 0.01 in Figs. [7] are limited to the region gpo > 1.4 where the contributions of low- 
energy eigenfunctions are negligible. Second, due to cancellations in the integral flBBl) we are 
not able to compute directly the probability current at large r (see, e.g.. Figs. [5] where the 
exact distributions are obtained at r < 45). This prevents us from normalizing directly the 
g^ = 0.01 distributions which fall off slowly at large r. Due to this feature we consider the 
normalization-independent quantity p(T)/pniax in Figs. O The exact data at g^ = 0.01 in 
Figs. E] are obtained by extrapolating the corresponding distributions with the conjectured 
behavior fl32|) . After extrapolation we normalize the distributions and compute the values 
of (r) and a^. 

We stress that the limitations listed above are not present at g"^ = 0.1; in this case the 
exact distributions p(r) can be computed directly at all values of gpQ. 
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